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We present an in-depth analysis of the experimental estimation of cross sections in Surface En- 
hanced Raman Scattering (SERS) by vibrational pumping. The paper highlights the advantages 
and disadvantages of the technique, pinpoints the main aspects and limitations, and provides the 
underlying physical concepts to interpret the experimental results. Examples for several commonly 
used SERS probes are given, and a discussion on future possible developments is also presented. 
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I. INTRODUCTION 

In the last decade Surface Enhanced Raman Scattering 
(SERS) 1] has made rapid progress towards applications. 
With a sensitivity rivaling fluorescence in some cases, and 
a much higher structural specificity, SERS is a highly at- 
tractive technique, being developed simultaneously with 
the field of plasmonics. Progress towards different uses 
of SERS in practical applications has been steady. New 
substrates including: arrays of inverted pyramids|2|, sil- 
ver pillar and torroid arrays [3, adaptive silver films 0, 
and metallic nano-shells[3|, have been demonstrated. In 
addition, many molecules relevant to a myriad of applica- 
tions such as gluco se Ig , p roteins H, DNAj^, a wide range 
of medicinal drugsj9[ll(lllll|. and substances for forensic 
science 12], have been characterized. On the other hand, 
the understanding of some fundamental aspects of the 
phenomenon is still incomplete, and in some cases con- 
troversial. Although a full understanding is sometimes 
not essential for the development of applications, there 
can be no denying the fact that better comprehension is 
desirable. 

Optical pumping of vibrational modes was first sug- 
gested in 1996 based on the observed dependence of 
the anti-Stokes/Stokes (aS/S) ratio with incident laser 
power J13J. and it was suggested that this could be a 
possible tool to estimate SERS cross sections. However, 
the interpretation of the experimental results have been 
the subject of considerable debate in the literature, with 
many authors simply denying its existence and attribut- 
ing the experimental observations to either laser-hea ting , 
resonance effects, or combinations thereof [Til llSl Ha. ItTJ . 
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Over a series of previous papers we have suggested the 
investigation of the aS/S-ratio as a function of tem- 
perature (r) as an alternative method for the obser- 
vation of pumping [iM, Il9l l2(]| . In a recent paper, we 
provided definitive evidence for the existence of vibra- 
tional pumping under SERS conditions thereby bringiirg 
to conclusion one aspect of this complex problem 20|. 
Here we address another aspect, namely the problems 
associated with estimating the SERS cross section us- 
ing this method, and we extend our analysis to include 
different analytes and laser excitations. We also high- 
light several outstanding issues which were not consid- 
ered previously J23. Obtaining good estimates of SERS 
cross sections is in general an extremely hard problem 
and has been a longstanding ambition of the SERS com- 
munity for reasons that go from the purely applied (quan- 
tification of signals) to the more fundamental (compar- 
isons of theoretical electromagnetic enhancement factors 
with experiment). Any method that can produce a stan- 
dard protocol for the estimation of SERS cross sections 
is, accordingly, of great interest and an effort to under- 
stand its principles and limitations is required. 

Before we go into the details, we briefly review the 
necessary concepts for understanding aS/S-ratios and vi- 
brational pumping in SERS. The method we are going 
to present has not been widely used and therefore we 
review the historical background and also its peculiar- 
ities and limitations. The discussion has some natural 
overlap with our previous papers [13, 13123 but it is pre- 
sented here for the sake of completeness, the convenience 
of the reader, and future reference for forthcoming work 
in progress. This will inevitably result in a somewhat 
lengthy introductory section where the principles of the 
method are laid down. The next sections are fully de- 
voted to this and are followed by a section at the end with 
a practical experimental demonstration of SERS cross 
section determination for several standard probes. 



II. A BRIEF DESCRIPTION OF SERS 
PUMPING 

A. The Anti-Stokes/Stokes ratio 

Let us consider the different contributions to the popu- 
lation of a single vibrational level at temperature T dur- 
ing a SERS measurement; we can identify two main con- 
tributions: (i) the laser itself which pumps vibrations 
through Raman processes with a rate proportional to 
its intensity (7^) and to the Raman-Stokes cross section 
(CT5), and (m) thermal excitation and relaxation. Vibra- 
tions remain in the level with a finite lifetime r, which 
encompasses all possible relaxation mechanisms -such as 
intramolecular vibrational relaxation IVR (anharmonic 
processes) or external relaxation mechanisms. There are 
other possible secondary mechanisms, such as relaxation 
through an anti-Stokes Raman processes, or excitation 
to higher vibrational levels. It is possible to write rate 
equations for the detailed dynamics of vibration popula- 
tions in such a system 'l(^. In the regime of weak pump- 
ing (which is the only case considered here), where the 
vibrational population remains small, n ^ 1, it is suffi- 
cient to consider only the two main mechanisms, and the 
rate equation for n can then be written as: 



Stokes intensities (p), we have: 



dt Tiujt, t r' 



(1) 



where 175 is the Raman-Stokes cross section, /^ the in- 
tensity (power per surface area) of the laser, hujh the 
energy of an exciting photon (n^, = Il/Tiujl is the num- 
ber of incident photon per unit time and surface area), 
and huju the energy of the vibration. The first term on 
the right is the number of vibrations per unit time being 
pumped into the level by the action of the laser, while the 
second and third terms are the contributions of thermal 
excitation and population relaxation, respectively. In the 
steady state dn/dt = 0. Generally, as is very small, so 
when II is small the pumping contribution is negligible 
and the vibrational population is dominated by thermal 
effects; i.e. n = exp{—hujy / ksT) is given by a Boltzmann 
factor. When II is increased (so that there is a signifi- 
cant pumping contribution to the population of the level) 
n becomes: 
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where the pumping term is clearly distinguished from the 
thermal contribution. 

In addition, for an ensemble of N molecules, the 
Stokes-Raman signal is given by Is — NasiL while the 
anti-Stokes signal is given by laS — nNoasIh, leading 
to: 



laS-i^-^+e-^-'^'^ANa^sh 
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where A is the asymmetry factor (extensively discussed 
in Refs. [H, [lE llal ) ■ ^ includes not only any possible 
difference between anti-Stokes and Stokes cross sections 
arising from resonance effects (due either to plasmon res- 
onances or to resonant Raman scattering for a resonant 
analyte), but also the standard wavelength dependence 
of Raman processes ^cj^ = {ijJl+uJi,Y /{ujl~uj^Y . p can 
be expressed as: 



p = A 



Tcrsh 



^-hui„/kBT 



hujj 



(5) 



Moreover, taking the ratio of the anti-Stokes to the 



This simple model shows that the Stokes-Raman sig- 
nal always remains linearly dependent on I^. The anti- 
Stokes signal also shows a linear dependence with power 
when pumping is negligible. But when pumping domi- 
nates over thermal effects, then las varies quadratically 
with power while the aS/S-ratio is linearly dependent on 

The original work by Kneipp et. al.[l3 used rho- 
damine 6G (RH6G) and crystal violet (CV) (two com- 
monly used SERS active dyes) under 830 nm excitation 
at room temperature (RT). The argument for pumping 
was primarily based on two observations: (i) the aS/S- 
ratios were shown to be larger than expected for a Boltz- 
mann factor, and (ii) the power dependence of the signal 
at RT (where the thermal contribution is likely to dom- 
inate, or be important) was quadratic. This observation 
of a non-linear dependence of the anti-Stokes intensities 
with power (resulting in a linear dependence in the aS/S- 
ratio) was thought to provide strong evidence for vibra- 
tional pumping in SERS. Further measurements on car- 
bon nanotubcs,.17j and a DNA-base|2j| supported these 
initial results. The arguments against this original inter- 
pretation are summarized in the following subsection. 



B. Resonance and heating effects 

Haslett et. al.jT^ were the first to seriously question 
the existence of vibrational pumping as revealed in the 
original studies. Extensive measurements were made un- 
der similar conditions with a number of both resonant 
and non-resonant molecules. They observed an anoma- 
lous ratio for all the resonant molecules tested which was 
independent of power until photo-bleaching occurred. No 
anomaly was observed in the case of the non-resonant 
molecules. It was concluded that the anomalous ra- 
tios observed in the original papers were the result of 
"hidden" resonances rather than pumping[23|. Brolo et. 
al. [13 also concluded that the observed anomalous ratio 
could be explained by resonances. These resonances are 
accounted for by the asymmetry factor A in the model 



of the previous section. In the absence of pumping, the 
aS/S-ratio is predicted to be: 



p — A exp 






(6) 



It is clear that the A factor can indeed result in anoma- 
lous ratios (different from the Boltzmann factor) , even in 
the absence of pumping. 

Furthermore, it was suggested that the observed power 
dependence was not the effect of pumping but rather laser 
heating. The effect of heating on p in the absence of 
pumping can be simply understood by including it in 
Eq. ©. The real temperature of the probed molecule, 
T/j, may be different to the nominal temperature T be- 
cause of heating effec ts^ either in the SERS substrate or 
in the molecule itself [la| ■ By writing T^ — T + AT and 
assuming that AT ^ T we can obtain the corresponding 
ratio, ph, by expanding Eq. |H|as: 



Ph = A exp 



huji, \ / huji, AT 

"k^J^'^^ykBT T 



(7) 



This expression first shows that heating can also result in 
an anomalous aS/S ratio, even when A « 1. Moreover, it 
can have an important impact on the power dependence 
of p, because AT should increase with II- We can as- 
sume in a first approximation that AT « alL , since heat 
diffusion and transfer are linear problems. The above ex- 
pression shows that laser heating should then result in 
an exponential increase of ph with I^. However, in many 
cases of interest, the argument in the second exponential 
is small compared to 1, and Eq. O can then be further 
expanded to give: 



Ph= A exp 
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A linear dependence of p with I^ (or equivalently a 
quadratic dependence of laS with I^) can then equally 
be the result of conventional heating effects or pumping. 
Such an observation is therefore insufficient to demon- 
strate the presence or not of vibrational pumping. We 
also note that heating effects can also explain the mode- 
dependent behavior observed in the original work [l3J, 
as the argument of the second exponential in Eq. ^ 
depends on the mode energy Tiuj^. 



C. Temperature dependence of the aS/S ratios 

Most studies have concentrated on finding evidence for 
SERS pumping at room temperature (RT), in particu- 
lar by studying the power dependence of the aS/S ratio. 
It is true that the huge enhancements in SERS condi- 
tions greatly increases the contribution from pumping, 
but may also contribute to an increased heating (in par- 
ticular directly in the probe molecules). At room tem- 
perature (RT) and above, the thermal contribution to 



p is relatively large and in many cases dominates, such 
that p w Aexp(— TiWiy/fc^T). The study of vibrational 
pumping at RT therefore involves measuring small de- 
partures from an already existing (large) thermal popu- 
lation. Moreover, as the previous discussion has shown 
it is also extremely difficult to distinguish the relative 
contributions of pumping and heating at RT. 

A more practical approach is to study the temperature 
dependence of the aS/S ratios, and in particular the low- 
temperature regime where thermal effects are expected 
to be completely negligible compared to any other mech- 
anisms, and in particular pumping. The temperature 
dependence can be directly studied using the model pre- 
sented so far. Figuren]summarizes the different scenarios 
for aS/S-ratios as a function of temperature. The solid 
lines show the variation of the ratio with T when pump- 
ing is absent; i.e. an exponential decrease as T decreases, 
following the Boltzmann factor exp{—hujy/kBT). The 
dashed lines show the case when pumping is present. At 
high T's the ratio is approximately the same as if pump- 
ing were absent. This is the thermally- dominated regime. 
As T is decreased, there is a cross-over to a pumping- 
dominated regime, where the ratio reaches a plateau. We 
refer to the cross-over temperature between these two 
regimes as Tcr] indicated by arrows in Fig. Q] The cross- 
over occurs when pumping and thermal terms in Eq. Q 
are comparable, which occurs for: 
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The cross-over occurs at a larger temperature for higher 
powers /i, for which pumping is stronger, but also for 
higher energy peaks (with a larger Tioj^), for which ther- 
mal excitation is weaker. 

In the pumping-dominated regime, the aS/S-ratio is 
constant and equal to p = AraslL/i^i^JL)- Il/{^l) 
can be estimated to a high degree of accuracy for a given 
experimental setup. If A = I (no asymmetry between 175 
and fJas) we can therefore deduce the product ras from 
the plateau in the aS/S-ratio below Tcr. In general, and 
in particular under SERS conditions, A ^ 1, but its value 
can be determined from the exponential dependence in 
the thermally-dominated regime (above Tcr). In practice, 
a fit of the experimental data over the whole temperature 
range with two parameters enables us to determine both 
A and ras [23 • Nevertheless, the problem of estimating 
T still needs to be circumvented to determine the SERS 
cross-section as itself. 

By reducing T, the contribution from pumping be- 
comes dominant, making measurements much easier, and 
more reliable. It also simplifies greatly their interpreta- 
tion. But the main advantage is that it enables one to 
unambiguously rule out heating effects as an alternative 
explanation. As discussed, the appearance of plateaus 
for T < Tcr in Fig. ^ is a clear characterization of the 
pumping-dominated regime. Note that it does not mean 
that heating is absent, but simply that its effect on aS/S- 
ratios is negligible (relative to the pumping contribution). 



Any heating effects may possibly affect the value of Tcr, 
but as long as a plateau is observed, vibrational pump- 
ing must occur. Moreover, the value of p in the plateau 
region (from which we will infer SERS cross-sections) is 
independent of heating. 

Studies of the temperature dependence of p therefore 
provide a more conclusive demonstration of SERS pump- 
ing. However, the technique is not without its drawbacks. 
In order to extract SERS cross-sections for the method 
discussed here there are several issues that need special 
attention as described below. 



III. PRACTICAL ESTIMATION OF 

CROSS-SECTIONS FROM VIBRATIONAL 

PUMPING 

Several issues need to be considered in the practical de- 
termination of SERS cross sections via vibrational pump- 
ing. Not all of them can be resolved satisfactorily and we 
have to include some ad-hoc approximations in order to 
obtain estimates of the cross sections. While there is 
nothing intrinsically wrong with this, it is necessary to 
be aware of the range of validity of the estimates in order 
to be able to comprehend the current limitations of the 
method as proposed. 



A. SERS cross-sections 

It is worth mentioning here that in most situations in 
Raman (and in SERS) we talk plainly about the cross 
section a while in reality there are at least three possi- 
ble cross sections we can refer to: (z) radiative, (ii) non- 
radiative, and (iii) total. In the case of vibrational pump- 
ing, for example, the rate at which vibrations are pumped 
into the molecule is proportional to the total cross sec- 
tion. But some of the re-emitted photons might not be 
observable in the far-field because they are absorbed by 
plasmon resonances in the vicinity of the probe. This is 
what makes the distinction between radiative and non- 
radiative cross sections. This could lead potentially to 
inconsistencies among cross sections for different modes 
if they happen at frequencies where the non-radiative 
contributions are different. We will not insist on this 
distinction in the following, but it is worth keeping in 
mind that as far as the cross sections is concerned what 
we measure is the "radiative" part of an effect (pump- 
ing) which depends on the "total" (radiative plus non- 
radiative) cross section. This distinction could under 
certain circumstances be important . 



B. Vibrational lifetimes 

We summarize in this subsection the different aspects 
of the problem of lifetime estimation for the purpose of 
obtaining SERS cross sections. 



• Ultimately, the observation of a cross-over from 
a thermally-dominated to a pumping-dominated 
aS/S-ratio leads to an estimation of both the asym- 
metry factor A and the product ras- An estima- 
tion of as itself requires the knowledge of r, as 
pointed out before. The latter cannot be directly 
obtained, in general, from the SERS spectrum po- 
tentially leading to a serious limitation of this tech- 
nique. In the first report of cross section estimation 
via vibrational pumpingjl^ a lifetime of r ~ 10 ps 
was simply estimated without any reference to the 
experimental data. 

In order to improve upon this we need to make a series 
of assumptions. The validity of these assumptions needs 
to be evaluated on a case-by-case basis and the results ob- 
tained for the cross section have to be interpreted within 
these approximations/assumptions. Let us first review 
the problems associated with trying to extract lifetimes 
from the SERS spectra. 

• Raman peaks in molecules have lineshapes which 
are seldom pure Lorentzians and have several con- 
tributions to their natural widths; a subject of a 
longstanding history in spectroscopy i235. The ob- 
served linewidth has typically contributions from: 
(i) groups of Raman modes piled up together in 
narrow energy regions; (ii) inhomogeneous broad- 
ening, (Hi) phase coherence relaxation, and (iv) 
population relaxation. The last two are also 
known as the off-diagonal (dephasing) and diag- 
onal (population) relaxation times in density ma- 
trix formalism 23] , respectively, and are part of the 
homogeneous broadening of the peak. These re- 
laxation times are different from mode to mode 
and cannot be easily separated from the observed 
linewidth in a plain SERS spectrum. The popula- 
tion relaxation is the one we need (t) for the esti- 
mation of as- 

• The presence of multiple modes contributing to a 
peak is a drawback that can be overcome in many 
situations. In the case of dyes it is not unusual 
for several Raman active modes, which are closely 
spaced in energy and have different cross-sections, 
to contribute to the observed peaks. This will natu- 
rally result in distortions of the lineshape and prob- 
lems with any estimate of the lifetime based on the 
width of the peak. This can be overcome by the 
use of well characterized vibrations which are rel- 
atively isolated from the others and even (if pos- 
sible) complementing the information by density 
functional theory (DFT) calculations of the Raman 
spectra|2^ ^^ to support the selection of those 
modes that are the best candidates for a reliable 
estimation of t. 

• If the broadening is homogeneous |2q . the only re- 
maining contributions to the linewidth are the diag- 
onal and off-diagonal relaxation times. Molecular 



vibrations have typically very strong anharmonic 
couplings with other vibrations in the molecule 
and contributions to inhomogeneous broadenings 
are mostly negligible. Even in SERS experiments 
where the single molecule limit is approached |27| 
(which would be more sensitive to inhomoge- 
neous broadening) the full width at half maximum 
(FWHM), r, of the peaks does not change by more 
than 1 — 1.5 cm^^ in peaks with a typical F of 
~ 15 — 20 cm"^. The validity of this assumption 
(whether inhomogeneous broadening is important 
or not) has to be assessed obviously in each specific 
situation. 

• The separation between diagonal and off-diagonal 
relaxation times cannot be achieved by a simple 
SERS spectrum. It is always possible to invoke a 
different type of spectroscopy (time resolved) [23 , 
but this may not be feasible in most situations. 
One simple possibility (as done in Ref. §31) is 
to estimate the population relaxation by using the 
FWHM (F) of the peaks via r ~ h/T] i.e. essen- 
tially ignoring the contribution of the dephasing 
(off-diagonal) relaxation to the width. This will 
produce an underestimation of r and an overesti- 
mation of as- In addition, this typically produces 
cross-sections for different peaks which are not en- 
tirely consistent with the relative (integrated) in- 
tensities among peaks; thus showing an intrinsic 
inconsistency in the estimation of the r's (which 
can be of the order of a factor of 2 to 10 for the 
examples we examined). It is generally difficult to 
make a simple estimation that will be valid for all 
modes. It is quite clear that for the cross sections 
to be meaningful (and if they are not affected by 
different amounts of non-radiative processes) the 
relative cross sections among peaks must be in ac- 
cordance with the relative integrated intensities ob- 
served in a normal SERS spectrum. 

• Population relaxation is caused mainly by the an- 
harmonic coupling to the "thermal bath" . This 
bath includes all other modes in the molecule as 
well as those in the solvent or substrate. Gener- 
ally speaking, higher energy modes have shorter 
lifetimes due to a greater number of possible de- 
cay pathways. Intramolecular anharmonic decays 
are energy conserving, meaning that high energy 
modes have more possibilities for a decay than low 
energy ones. Population relaxation is also the only 
dynamical process that contributes to the vibra- 
tional linewidth in the limit T ^ 0. 

• It is important to realize that, in principle, if we 
find a single peak in the spectrum where the pop- 
ulation relaxation can he gained directly from the 
FWHM via r ^ h/T, then all of the cross sections 
of the other modes follow immediately through the 
relative integrated intensities of the peaks. When 



lacking a time resolved experiment to isolate the 
contributions properly, the problem comes down to 
a judicious choice of a mode where the width is 
dominated by relaxation population. A compromise 
then is to extract the values of the cross sections 
by means of the width of the highest (isolated) Ra- 
man active mode. A rule of thumb is that the lower 
the temperature and the higher the energy of the 
mode the more the linewidth will be dominated by 
population relaxation. We can then use one mode 
to estimate its as, and all the others follow auto- 
matically from this assumption. The cross sections 
are now consistent by construction and the method 
provides, in fact, the means to estimate the relative 
lifetimes of the modes. We call this procedure (i.e. 
estimating the r of the Raman active mode with the 
largest Raman shift, and making the other cross- 
sections consistent with this value through the rel- 
ative integrated intensities) the corrected lifetime 
method (CLM). 

• If we are comparing different substrates with the 
same analyte we could also directly compare the 
values of ras without any further assumption. If 
there is no reason to believe that relaxation popu- 
lation is different between the two substrates, the 
ratio of ras for two different .substrates provides a 
direct comparison of SERS cross sections. 



C. The asymmetry factor 

• Problems with measuring the asymmetry factor 
"A" can be as important as those to obtain a reli- 
able lifetime. Heating in the thermally-dominated 
regime can be a major problem. As demonstrated 
in Sec. II B, heating can be in part a source for a 
larger than normal A, if it is not properly identified. 

• Ideally, the best estimation of A can come from 
measurements at RT with long integration times 
and very low powers. The obtained A's should be 
more or less the same for a given analyte in the 
presence of the same metal (they are mainly re- 
lated to the local interaction of the dye with the 
metal_22j). 

• The product Aras is measured with high accuracy 
in the pumping dominated "plateau" at low tem- 
peratures. Accordingly, the best use of the tech- 
nique is to use one well-characterized analyte (with- 
out photo-bleaching if possible) as a means to com- 
pare the relative SERS cross-sections of the dye on 
different substrates of the same metal, by simple 
comparison of the plateaus. The product Aras is 
perhaps the best "objective" comparison between 
the cross sections of two substrates without any 
assumption. However, it is also important to note 
that the cross sections we are comparing are not the 



average cross sections but rather a biased average 
towards the highest cross sections in the sample, as 
explained in the next subsection. 



D. Which cross section do we really measure? 

There are additional complications with the estimation 
of the cross section. It is generally accepted that SERS 
signals are dominated by the presence of "hot-spots" or 
places with high local enhancements. The formulae pre- 
sented above made implicitly the assumption that the 
cross section is the same for all N molecules in the sam- 
ple. The question we want to address now is: do we 
obtain an estimate of the average cross section with this 
method? We shall show in what follows that what we ac- 
tually obtain from the experiment is an estimate which 
is heavily biased towards the sites with the highest en- 
hancements; i.e. the method provides an estimate of the 
enhancements in hot-spots. This is of course an advan- 
tage and a disadvantage at the same time. We show this 
explicitly in this section. 

We consider a number of molecules (z from 1 to N) 
in the sample. Each molecule can have a different anti- 
Stokes, {(Jas) ^^'^ Stokes (ag) SERS cross sections (ac- 
cording to the asymmetry factor A* = '^as/'^S ^'^^ ^'^^^ 
molecule). We assume for simplicity that the incident 
power /^ is the same for all. The rate equation for the 
average phonon population n' of a vibrational mode of 
one molecule, i, and its stationary state solution are given 
by expressions similar to Eqs. ^ and (|2Jl: Following 
the other formulae in the previous sections, the Stokes 
signal for this molecule is simply Ig = o-gl^, while the 
anti-Stokes signal is given by: /*g = n^al^glL- The anti- 
Stokes/Stokes ratio for this molecule, p*, is given by an 
expression similar to Eq. (|SJ). If all the molecules expe- 
rienced the same enhancements, then the total Stokes 
and Anti-stokes intensities would be Is = Nig and 
laS = Nllg, and the ratio would be p = pi. However, 
if the molecules have different cross sections (SERS en- 
hancements) , which is a much more realistic assumption, 
we then have: 



semble averages. The aS/S-ratio is then: 
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All the averages here (...) are the usual statistical en- 



P = las/h 









(13) 

Comparing this expression with that obtained previ- 
ously in Eq. Q, what we called A is now replaced by: 



A' 






(14) 



Note that this is not exactly the average of the A"s, i.e. 
{A). Moreover, what we measure (instead oiras) is now: 



TCTc 



{<^as) 



(15) 



To understand the meaning of ag {E for Ensemble), we 
need to understand the sources of non-uniformity. The 
main variation is in the cross-sections, because of the 
large differences in SERS enhancements at different lo- 
cations on the SERS substrate. We typically have val- 
ues which arc, say, 10"^ — 10^ larger at hot-spots than 
in other places. However, for a given molecule, the ra- 
tio A' = crl^s/o-g should not vary as much (at least not 
by more than one order of magnitude). In fact, this ratio 
can in a first approximation be taken as a constant, which 
is an intrinsic property of the adsorbed molecule/metal 
complex [2^. Assuming that A^ = A is the same for all 
molecules, we simply get that Ae = A, i.e. the A we 
measure is the correct one. Moreover, using cr^g = ^'^s^ 
we can express a^ as a function of (Tg only as: 
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If the distribution of erg's were fairly uniform (for ex- 
ample a Gaussian around an average value with a small 
standard deviation) then a^ would be a good estimate of 
the average of the distribution (slightly overestimated). 
But this is far from reality in SERS conditions. A more 
realistic situation is to have a small number of molecules 
at a hot-spot (HS), and a large number of molecules 
at non-HS positions. For the sake of argument, let us 
demonstrate the effect on this average by considering 10'^ 
molecules with one molecule at a HS and 999 molecules 
at places with much lower enhancements. The expression 
given above would give: 






1/1000 (l(afS) 



999((7NHS)2) 



1/1000 (l((7HS) + 999(fTNHS)) 



(17) 



If there is a difference of a few orders of magnitude be- 
tween the cross sections cr^ ^ and cr^ ^^ , the numerator is 
then easily dominated by the one molecule at the HS (due 
to the square). For the numerator, it is not as clear-cut, 
it all depends on how much stronger the HS is and how 
many molecules are at non-HS positions. But overall, the 



denominator is likely to be dominated by Cg, possibly 
slightly larger if the second contribution is not negligible. 



(Jg , or perhaps slightly smaller. 



This shows that ag 
A more quantitative argument is only possible if we had 
a more realistic distribution of SERS enhancements for 
a given substrate, but the qualitative conclusion is that 
pumping experiments provide a good lower estimate of 
the cross-section of the few molecules experiencing the 
highest enhancements. 



E. Photo-bleaching 

The last aspect we want to briefly touch upon in 
this analysis of cross section estimations via vibrational 
pumping is photo-bleaching. The exact mechanisms of 
photo-bleaching under SERS conditions are still poorly 
understood and would deserve a full study in itself. 
Photo-bleaching is particularly important in the method 
we are describing here because molecules sited at HS's 
(i.e. those exposed to the largest enhancements) are 
affected to a greater degree. At high power densities, 
for example, the population of molecules at HS's will 
be bleached at a faster rate than the average and this 
will change the measured cross section. Studying p un- 
der photo-bleaching conditions could, in principle, tell 
us something about this distribution. There is robust 
experimental evidence that photo-bleaching can compli- 
cate the analysis, produce experimental artifacts, and -to 
a large extent- even decide the answer we obtain from 
a specific experiment. This is added to several exper- 
imental complications in aS/S-ratios which include the 
fact that sometimes (depending on the dispersion of the 
spectrometer) the anti-Stokes and Stokes sides cannot be 
measured simultaneously. This implies a delay between 
the two measurements and therefore a difference in expo- 
sure times to the laser. By choosing a certain power level 
we are, in a way, selecting the population of HS's that we 
are going to measure. We shall come back to this prob- 
lem in the discussion of the results. A rule of thumb is 
that larger laser spots with low power densities (but suf- 
ficient to produce detectable pumping), short integration 
times, and non-resonant (with the dye) laser excitation, 
are in general preferable for a more reliable estimate of 
SERS cross sections. All possible measures should be 
taken (including the type of scanning used) to address 
and minimize the undesirable effects of photo-bleaching. 



IV. EXPERIMENTAL 

We turn now to an experimental demonstration of the 
principles underlined above. SERS measurements have 
been performed on substrates formed by dried Ag col- 
loids on silicon. The colloids were prepared using the 
standard Lee and Meisel technique 28]. SERS active 
samples were prepared by mixing the colloids with a 
20 mM KCl solution in equal amounts. The analytes 



were then added to give a concentration of 10~^M in 
each case. A small amount of this solution was then 
dried on to a silicon substrate. The investigated ana- 
lytes are rhodamine 6G (RH6G), crystal violet (CV), and 
3,3'-diethyloxadicarbocyanine (DODC). Figure |2] shows 
the basic SERS spectra of the analytes. We shall show 
the explicit temperature dependencies of the aS/S-ratios 
for only a few modes (labeled in Fig. |2Jl while the ta- 
bles show additional data for other modes. Samples were 
mounted in a closed-cycle He-cryostat (CTI-Cryogenics) 
with temperature control in the range 10 — 300 K. Ra- 
man measurements were performed using several laser 
lines of a Kr"*" and Ar"'"-ion lasers which were focused 
to a 20 /im diameter spot. The signal was collected us- 
ing a high-numerical aperture photographic zoom lens 
(Canon, xlO magnification) onto the entrance slit of a 
high-dispersion double-additive UIOOO Jobin-Yvon spec- 
trometer coupled to a liquid N2-cooled CCD detector. 
Peaks were analyzed using standard Voigt functions with 
subtracted backgrounds. The variation in the aS/S-ratio 
with T for each mode was then fitted to Eq. [S] For this 
purpose, it is convenient to modify the expression to: 



ln(p) 



In 



-huj^/kisT 



(18) 



where a = ln(A), and b — tctsIl/TT'^l which are the 
(dimensionless) fitted parameters. The t's of the peaks 
(necessary to obtain the cross sections) are estimated in 
two different ways following the prescriptions in Sec. Ill 
B. 

As discussed earlier, it is of critical importance that the 
sample remains stable over the entire length of the exper- 
iment. A range of integration times were used depend- 
ing on the specific sample and measurement to achieve a 
good signal to noise ratio for different peaks, but power 
densities were kept to a minimum compatible with the 
observation of pumping and exposure to the beam was 
minimized as much as possible in between changes of T . 
We used ~ 30 — 50 mW spread over the 20 /im diameter 
spot. As the rate of photo-bleaching may change with 
T 29] a suitable power density to ensure reasonable sta- 
bility must be decided at RT, where the effect is greatest. 
The relatively large spot and low power densities have 
the triple advantage of (z) reducing photo-bleaching to a 
negligible level, {ii) reducing any indirect laser heating 
effects and {Hi) improving the averaging over cluster ge- 
ometries. Still, even under low power density excitation, 
the choice of laser line can have an important effect on 
the measurement. An explicit example is shown in Fig. |3| 
for RH6G under 514 nm excitation. As explained in the 
caption, this observation is compatible with a bleaching 
of the HS's with largest cross sections. Laser lines that 
are strongly resonant with the probe molecule 22] should 
be avoided so that the photo-stability of the analyte and 
the accumulated exposure to the laser do not become 
an issue in the interpretation and analysis of the data. 
For that reason, we only concentrate hereafter on results 
obtained with the 647 and 676 nm lines of a Kr+-laser. 
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These two lines are close enough to the visible (to profit 
fully from SERS enhancements) but, by the same token, 
they do not produce excessive photo-bleaching for the 
dyes under consideration here. 



V. RESULTS AND DISCUSSION 

Having all the theoretical tools from the previous sec- 
tions, we can now easily scan through several experimen- 
tal results that demonstrate the method in practice. Fig- 
ure^ shows the aS/S-ratios (In(yo)) vs. T, as measured 
using the 676 nm line, for four modes of each of the inves- 
tigated analytes. The solid lines represent the best fit to 
the experimental values using Eq. ^| A series of small 
imperfections can be seen in the data for some modes, 
but overall the behavior is very well represented by Eq. 
1181 with RH6G being the best example. The arrows in- 
dicate the cross-over points (Tcr) between the thermally- 
dominated and pumping-dominated regimes for all the 
investigated modes. Tcr occurs at higher T's for higher 
energy modes, as expected. The fact that pumping is ob- 
served for each of the modes indicates that the effect is 
fairly general and should be observable for a large num- 
ber of analytes for as long as 175 is high enough to make 
pumping observable at moderate power densities |8n|. 

In order to transform these data into estimates of the 
cross sections we need to address the different meth- 
ods of estimating r. Table HI shows the experimental 
results for RH6G taken with the 676 nm line. For the 
sake of comparison, we display the cross section by using 
two different estimations of the r's. This gives a better 
idea of how reliable these numbers are and how differ- 
ences among them should be interpreted. Table Q] shows 
both: (i) the extraction of r directly from the widths via 
T ~ h/T (this normally renders cross sections which are 
not fully consistent with the relative integrated intensi- 
ties of the peaks), and (m) the extraction of r from the 
highest Raman active mode and the subsequent adjust- 
ment of the other cross-sections to agree with the rela- 
tive integrated intensities (as described in Sec. Ill B). A 
column with relative intensities and cross-sections among 
modes is also provided for completeness. The anti-Stokes 
cross sections are obtained through the Stokes ones via 
the fitted value of the asymmetry parameter A from Eq. 
[THl These results prove the point that the estimation 
of the lifetime is a crucial step in this method to trans- 
form the experimental values into an estimation of the 
cross sections themselves. The values should be inter- 
preted within this assumption, and their validity should 
be assessed on a case-by-case basis. The preferred (more 
consistent) method is of course the one that ensures the 
relative cross sections agree with the relative intensities 
of peaks (CLM in Sec. Ill B). 

Finally, we can discuss the effect of estimating the 
cross sections for two different excitation wavelengths. 
The experiment is somewhat restricted in the choice of 
excitations (to be compatible with photo-stability) but 



still provides a hint of resonance contributions to the 
cross sections with some limitations. Table UTI shows the 
cross-sections obtained for RH6G, DODC, and CV for 
647 and 676 nm excitations for different modes. All the 
cross sections were obtained by the CLM analysis. It 
is worth noting that in this latter method we only need 
(in principle) one temperature dependent aS/S-ratio for 
the highest Raman active mode (judged to be dominated 
by population relaxation). This together with the plain 
SERS spectrum of the dye is enough to estimate all the 
cross sections of the modes via the relative integrated in- 
tensities. The results will be consistent (by construction) 
but not necessarily accurate (depending on the one mode 
that has been chosen to estimate the initial cross section) . 
The main conclusions of these experiments can be sum- 
marized as follows: (i) RII6G shows larger cross-sections 
at 647 than at 676 nm. This is consistent with a resonant 
increase known to exist in RII6G towards the yellow part 
of the spectrum 22). (ii) On the contrary, both DODC 
and CV experience larger cross-sections at 676 nm, con- 
sistent with the maximum of the absorption of these dyes 
being at longer wavelengths with respect to RII6G and, 
therefore, having a much more red-shifted resonant in- 
teraction than RH6G|2a|. {Hi) Although there is some 
variations from mode to mode, it is generally true that 
CV experiences the largest cross-sections in this near-IR 
region. 

The numbers in tables ^ and ^ deserve a special com- 
ment regarding their interpretation. For a start, the main 
source of error in these numbers (as pointed out in the 
caption of table P) is systematic rather than statistical, 
and comes from the different methods used in the esti- 
mation of T. Error bars for these systematic errors are 
not easy to obtain from the data, for they depend on 
the relative contribution to the FWHM of population re- 
laxation and dephasing, which is unknown and changes 
from mode to mode. The number of significant figures is 
decided then by the range in which these numbers span 
for a given group of modes and laser excitation. It only 
makes sense to compare statistical errors (from the prop- 
agation of errors in the parameters of the fits and the 
experimental magnitudes) for two cross sections with the 
same systematic errors. This is the situation, for exam- 
ple, of the comparison between the method presented in 
this paper and the power-dependence method (quadratic 
dependence of the las signal as a function of power) as 
reported in Ref . 20] . In this latter case the same system- 
atic error coming from the estimation of r affects both 
results and, hence, the statistical error bars are meaning- 
ful to claim consistency or inconsistency between the two 
values. In any other case, the overriding source of error 
in these estimates of the cross sections is systematic. 

The order of magnitude of these numbers deserves a 
special comment too. According to the results in Sec. 
Ill D, we expect these numbers to represent something 
close to the maximum available enhancement in these 
samples. In addition, the fact that r could be in general 
underestimated by any of the methods presented here 



leads automatically to an overestimation of a. Last but 
not least, the values of OaS f-re obtained from those of 
Us via the asymmetry factor A. As a consequence of all 
these, we expect the values to be: (i) At the higher-end 
of what is known for HS's in SERS for ag and, (ii) with 
a larger error for the anti-Stokes values which use the 
estimation of the asymmetry parameter (^) two times; 
once to obtain ag itself through the fit of ln(p) vs. T, 
and twice to obtain tJaS from as- Indeed, Table HTl shows 
that cross-sections are at the higher-end of what is com- 
monly accepted for SERS enhancements, a^s is bound 
to be less exact than as because it depends crucially 
on the asymmetry parameter and its associated errors 
(treated in Sec. Ill C). The most consistent set of val- 
ues are obtained for as using the CLM-analysis (which 
should be the preferred option) and are all of the order 
of ^ 10^^^ cm^. In exceptional cases, cross sections of 
the order of ~ 10~^^cm^ are obtained, which are still 
expected if we take into account that the estimation is 
for HS's and that the values are possibly overestimated 
by the choice of t. Unrealistically large values of as be- 
yond that should be interpreted as a limitation of the 
technique, in general, and should be tracked down to a 
problem in the assignment of r for that mode. Even for 
these cases, it is still possible to compare Tas (or Aras^ 
for the same modes in different substrates providing a 
relative comparison of cross section performances at HS's 
without any additional assumptions. This comparison is 
independent of the number of molecules involved in the 
signal, due to the self-normalizing nature of the method. 
Table ITTI summarizes, therefore, a typical example where 
the performance of different dyes can be tested through 
appropriate estimates of the cross sections derived from 
the method presented in this paper. 



VI. CONCLUSIONS 

We have presented an in-depth discussion of the state- 
of-the-art of SERS cross section estimation via vibra- 



tional pumping and have given several experimental ex- 
amples of it. Despite the approximations needed and the 
intrinsic problems of the technique, it is still an excellent 
tool (and in many cases possibly the only tool) to extract 
an estimation of the SERS cross-sections in situations 
where the number of molecules and hot-spots present in 
the sample are not known or are difficult to estimate. 
The self-normalizing nature of the method (with respect 
to the number of molecules involved) is certainly a major 
advantage. In addition, we feel that these latter experi- 
mental and theoretical developments have moved forward 
to establish vibrational pumping as a real phenomenon in 
SERS after almost a decade of controversy about its very 
existence. One outstanding issue that surely deserves fur- 
ther investigation is the effect of photo-bleaching under 
SERS conditions (in particular in HS). As we showed in 
Sec. HI D, the fact that this technique produces an esti- 
mate of the cross sections which is heavily biased towards 
HS's provides a unique opportunity to explore their char- 
acteristics. This is particularly true in the comparison 
among different substrates (in which the photo-stability 
of the analyte plays a special role) . A forthcoming study 
on the topic is in preparation and will be published else- 
where. 
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FIG. 1: Schematic diagram showing the temperature depen- 
dence of the anti-Stokes/Stokes ratio (p) for several different 
Raman modes. We have assumed here A = 1 for simplicity 
and have used three characteristic Raman modes of RH6G for 
the example. The solid lines represent the ratio when there 
is no pumping contribution (Boltzmann factor). The dashed 
lines show the result of including pumping as given by Eq. 
1^ There are two regimes, one at high temperature where 
the thermal contributions dominate and the ratio is similar 
to what is expected if no pumping were present. The second 
occurs at low T's where the contribution of pumping becomes 
significant relative to the thermal contribution. We refer to 
the cross-over point between these two regimes as the Tcr (in- 
dicated by the arrows in this figure). Tcr generally occurs at 
higher T's for higher energy modes due to the different influ- 
ence of the exponential factor for different mode energies in 
Eq. 13 See text for further details. 
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FIG. 2: SERS spectra of the three analytes used in this pa- 
per taken with a 633 nm laser: crystal violet (CV, top), 3,3'- 
diethyloxadicarbocyanine (DODC, center), and rhodamine 
6G (RH6G, bottom). The temperature dependence of the 
aS/S-ratios of the peaks labeled with "•" are explicitly shown 
in Fig. 0] Table UTI includes, in addition, estimations of cross 
sections for additional peaks not labeled here. 
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TABLE I: Values of 6 = raslL/JiLOL, t, relative integrated intensities, and cross-sections of the Stokes modes. The values of as 
and aaS are shown for two analysis schemes - (i) r's obtained from the FWHM of the peaks (ignoring dephasing contributions), 
and {a) corrected r's to account for the relative integrated intensities. All data are for RH6G using the 676 nm laser line 
as excitation. The "Rel. Is" column is obtained directly from the relative integrated intensities of the peaks in the SERS 
spectrum. Cross sections for the different estimates of r are rounded to one significant figure decided by the range of variation 
of a over the modes. We do not quote error bars because the overriding cause of error in these numbers is systematic (rather 
than statistical) and comes from the estimation of r itself. The systematic errors are not easily obtainable from the data, for 
the exact amount of dephasing contribution to the broadening in each peak is unknown. 
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FIG. 3: Effect of photo-bleaching on the temperature depen- 
dence of the aS/S-ratio at 514 nm excitation. The measure- 
ment in (a) is performed at low power density (25 mW on 
a 20 ^m spot in diameter). The sample is then exposed to 
400 mW for 3min on the same spot and the measurement is 
repeated afterwards with the same original low power density. 
The decrease in the plateau of ln(p) to smaller values (and the 
associated shift of T^r to lower temperatures) in (b) is consis- 
tent with a bleaching of the molecules on the HS's with the 
highest enhancement in the first measurement. The dashed 
line in (b) is the value at which ln(p) tails in (a). By choos- 
ing a certain power level, we are selecting the population of 
HS's that will survive through the experiment. We are there- 
fore measuring a convoluted property of the substrate and the 
photo-stability of the probe. 
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FIG. 4: Anti-Stokes/Stokes ratios (ln(p)) as a function of tem- 
perature for four modes of each of the analytes investigated 
here (all taken using 676 nm excitation) . The arrows indicate 
the crossover points between the thermally-dominated and 
pumping-dominated regimes. Note that the cross over occurs 
at different temperatures for different modes, being at higher 
T's for larger vibrational energies. The solid lines represent 
the best fit to the experimental data using Eq. 1181 Table UTI 
presents an analysis of these data using different assumptions 
for the lifetimes. 
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124 


1.1 


3.9 


6.82 


1175 


19.5 


30.7 


18.0 


4.3 


31.8 


5.62 


1292 


3.7 


6.3 


8.24 


1.6 


4.8 


3.94 


1620 


6.4 


14.4 


9.99 


2.0 


32.7 


2.26 



TABLE II: Cross-section for the Stokes and anti-Stokes Raman modes of CV, RH6G and DODC extracted from b = raslL/hcuL 
obtained through fitting the experimental data to Eq. ll8l for both the 647 and 676 nm lasers. Only the lifetime of the mode with 
the largest Raman shift is estimated from its width and the other cross sections are made consistent with this determination 
via the relative integrated intensities (CLM- or "corrected r's" method in Sec. Ill B). See the text for further details. 



